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Abstract 

We study the relation between the causahty and positivity of energy bounds for Gauss- 
Bonnet gravity in an AdS^ background. Requiring the group velocity of metastable states 
to be bounded by the speed of light places a bound on the value of Gauss-Bonnet coupling. 
To find the positivity of energy constraints we compute the parameters which determine the 
angular distribution of the energy fiux in terms of three independent coefficients specifying 
the three-point function of the stress-energy tensor. We then relate the latter to the 
Weyl anomaly of the six-dimensional CFT and compute the anomaly holographically. The 
resulting upper bound on the Gauss-Bonnet coupling coincides with that from causality 
and results in a new bound on viscosity /entropy ratio. 
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1. Introduction and summary 

Recently an intriguing connection between positivity of energy flux and causality has been 
addressed in the context of higher derivative gravity [1-5]. Consider Gauss-Bonnet gravity 
with negative cosmological constant in five dimensions. Since AdS space is a solution of 
the equations of motion, one can hypothesize the existence of a dual CFT. The theory 
furthermore possesses exact black hole solutions which are asymptotically AdS. Brigante, 
Liu, Myers, Shenker and Yaida [2-3] considered gravitons propagating in these backgrounds 
and found long lived excitations which correspond to metastable states in the boundary 
theory. Restricting the group velocity of these states to be bounded by the speed of 
light (requiring causality of the boundary theory at finite temperature) places non-trivial 
constraints on the value of Gauss-Bonnet coupling A [2-3] . 

At first sight, a completely unrelated set of constraints was proposed for CFTs by 
Hofman and Maldacena [1]. By requiring positivity of the energy fiux measured in a 
collision of certain CFT states, they deduced a set of constraints on the quantities t2 and 
t4 which determine the angular distribution of the energy flux. The values of t2 and t4 are 
completely determined by the two- and three-point functions of the stress energy tensor. 
The constraints can be reformulated as bounds on the ratio a^^^^^ / c^^^^^ of the coefficients 
which appear in the Weyl anomaly of a four-dimensional CFT. In the supersymmetric case 
t4 = and the bounds are 

i < < ^ (1.1) 

These bounds have been verified for a number of interacting superconformal theories [6]. 
(More stringent bounds exist for jV = 2 superconformal four dimensional field theories 
[7].) In the non-supersymmetric case a^^^^^ and c^''^^) do not completely determine t2 
and t4, and therefore the bounds on a^'^^^^/c^''^^) are slightly relaxed. In [1] it has been 
noted that for Gauss-Bonnet gravity the lower bound in (1.1) precisely translates into the 
upper bound on A coming from causality^. Recently the upper bound in (1.1) has been 
reproduced by examining a differently polarized graviton (shear channel) [4] (see also [5]). 

At first sight this relation is puzzling. Why would the J\f = 1 superconformal result 
have anything to do with Gauss-Bonnet gravity? It is worth to note however that the 

^ The upper bound on A implies a new bound on the viscosity /entropy ratio in higher derivative 
gravity [2-3]. The original KSS bound r]/s > 1/4-ir [8] has been shown to be violated in a controlled 
setting [9]. Other recent work on the higher derivative corrections to 77/s includes [10-23]. 
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terms are absent in supersymmetric string theories, so one may fantasize that some 
supersymmetric string compactification would yield Gauss-Bonnet gravity as the low en- 
ergy theory. In this paper we investigate the state of the correspondence between causality 
and positivity of the energy flux in the case of Gauss-Bonnet (GB) gravity in an AdS^ 
background possibly dual to some six-dimensional CFT. An important feature of the six- 
dimensional CFT is that, unlike the four-dimensional case, the knowledge of the Weyl 
anomaly completely determines the values of t2 and ^4. In particular, we find that ^4 = 0. 
Note that vanishing t4 is a necessary feature of any supersymmetric CFT. We compute the 
upper bound on A by requiring both causality of the boundary theory at finite temperature 
and positivity of the energy flux. The two bounds coincide exactly and result in a new 
bound on rj/s. We also compute the lower bound on A. 

The rest of the paper is organized as follows. In the next Section we review Gauss- 
Bonnet gravity and asymptotically AdS black holes solutions. In Section 3 we study small 
fluctuations around these black holes, and flnd metastable states which can lead to causality 
violation. By requiring the group velocity of these states to be bounded by the speed of 
light, we flnd an upper bound on the value of A. In Section 4 we generalize the results 
of [1] to six dimensions. We compute the values of t2 and in terms of the parameters 
a, b, c which determine the two- and three-point functions of the stress-energy tensor. In 
Section 5 we compute the coefl&cients a, 6, c (and hence t2 and 14) in GB gravity. We do it 
by relating the values of a, 6, c to the coefficients of the B-type terms in the Weyl anomaly 
and computing the anomaly holographically. We then find the bounds on A and observe 
that the upper bound is precisely the same as the one found in Section 3. We discuss our 
results in Section 6. Some technical results appear in the Appendix. 

2. Gauss-Bonnet gravity 

Among gravity theories which involve higher derivative terms of the Riemann tensor in 
their actions, there is a special class which shares many of the properties of Einstein- 
Hilbert gravity. It is usually referred to as Lovelock gravity [24-26] and is the most general 
theory of gravitation whose equations of motion contain at most second order derivatives 
of the metric. Recently, the Palatini and metric formulations of Lovelock gravity have been 
shown to be equivalent [27] . 
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The Lovelock action for a d + 1-dimensional spacetime can be written as 




Here Ip is Planck's length, [|] denotes the integral part of |, Ap is the p-th order Lovelock 
coefficient and Cp is the Euler density of a 2p-dimensional manifold. In c? + 1 dimensions 
all Cp terms with p > [|] are either total derivatives or vanish identically. 

The Gauss-Bonnet action is the simplest example of a Lovelock action, with only the 
4-dimensional Euler density included. In the following we set d = 6 




Note that in eq. (2.2) we set Ip = 1, included a cosmological constant term A = — ^ and 
rescaled the Lovelock parameter by L^. The Gauss-Bonnet term £(2) is 

£(2) = RmnpqR^^^^ — 4:RmnR^^ + R^ (2.3) 

Equations of motion derived from (2.2) can be expressed in the following way 

1 1 

— -^Qmn^ + Rmn + Q^'Hj^j^^ = (2.4) 

with T-L^MN defined as 

^MAT — RmlpqRn ~ "^RmpRn ^ ~ '^RmpnqR^^ + RRmn (2.5) 
Eq (2.4) admits a solution of the form 

ds^ = -a'f{r)dt' + J^^+ J2d^lk (2.6) 

with dTi\ ^ the metric of a 5-dimensional manifold of constant curvature equal to 20/c [28- 
29]. Note that a in eq. (2.6) is an arbitrary constant which allows one to fix the speed of 
light of the boundary theory to unity. Given that we are interested in black hole solutions 
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with flat horizon, we set A; = in the foUowing. In this case the solution, known to be 
thermodynamicaUy stable, reduces to^ 



'l-4Afl-% 



(2.7) 



a = 



1^2 

lim -^f{r) 



1 + Vl - 4A 



1 — >oo r 

The horizon is located at r = r+ whereas the Hawking temperature of the black hole is 

,3 r_|_ 



T = a 



2L2 



(2.8) 



In the limit >0 one recovers AdS'j space from (2.7). The curvature scale of the AdS 
space is related to the cosmological constant via 



LacLS = ctL 



(2.9) 



3. Fluctuations 



In this section we study fluctuations around the black hole solution (2.7). In particular, 
we consider small perturbations of the metric /imat in the scalar channel (j) = hi2- In this 
case, the form of the perturbed metric is 
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ds' 



-a''f{r)dt' + 



dr' 



+ 



f{r) 



dx^ + 2(f){t, r, x^)dxidx2 



.i=l 



(3.1) 



Note that (j) only depends on the (t, r, x^) directions of spacetime and its Fourier transform 
can be written as 

0(t, r, 0:5) = j |^^(r)e--*+^«-^ k = (a;, 0, 0, 0, 0, q) (3.2) 

The equations of motions for (f can be found by substituting the ansatz (3.1) into (2.4) 
and expanding to linear order in the fluctuating field. The result is 



T2^"(r) + Ti(^'(r) + ro(^(r) = 



(3.3) 



^ Gauss-Bonnet gravity admits another AdS solution with = | [l — Vl — 4A] . This 
however unstable and contains ghosts [28]. 



IS 
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with 

To = 3ruj^ [-2r^ + XL^r'^f'{r) + 2rXL^f{r)\ + 

+ fa^L^fir) [-A\L^rf'{r) + 6r^ - \L^r'^f"{r) - 2\L^f{r)\ 
Ti = Za^L^rf{r) [r^f'{r) (-2r + \L^f'{r)) + QXL^f{rf (3.4) 

+r/(r) {8XL^f'{r) - lOr^ + rXL^f"{r))] 
T2 = Sa^L'^rfirf [-2r^ + XL\^f'{r) + 2rXL^f{r)] 

where primes indicate differentiation with respect to the variable r and we defined u = ujLP' 
and q = qL^. 

It is convenient to place this equation in Schrodinger form. To do this, we follow two 
steps: We first define a new function $(r) through 



ln^ = ln^+^/| 



(3.5) 



which brings eq. (3.3) into the standard form 
$"(r) + W{r)^(r) = W{r) 



To 


1 


'T[T2-T!,T^' 




T2 


~ 2 


[ n \ 





(3.6) 



Eq. (3.4) then allows us to express To(r) as To(r) = bi^{r)oo'^ + bq{r)q^ where 6a;(r), bq{r) 
are functions of the radial coordinate r alone. This implies that W{r) can be written as 



W{r) = ^^u' + + h{r) h{r) = 



1 


'T[T2-T!,T^' 




m 


2 


[ n \ 







(3.7) 



We then substitute $(r) with '^[r) defined as '^{r) = ^ ^{f) and subsequently make 
a coordinate transformation from r to y according to 



9ry{r) = \ 7^ 

V -^2 

Eq. (2.4) is finally expressed as 

-dl^+[q'cl{y) + Vi{y)]^^u^^ 

where 

o bn 



(3.8) 



(3.9) 



c„ = 



vfe) = -|%)-(|) s. 



T2 '[t2 



(3.10) 
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We are now ready to study the full graviton wave function (3.9) . Note that y{r) is a 
monotonically increasing function of r with j/— >0 at the boundary r >> r+ and j/— > — oo 
at the horizon r = r+. Vi{y) blows up as y~'^ for y^O. 

Following [2,3] we consider (3.9) in the limit g— )-oo. In this case, q^c^g{y) provides the 
dominant contribution to the potential except for a small region y > —k. It is therefore 
reasonable to approximate the potential with c^{y) for all y < and replace it with an 
infinite wall at j/ = 0. Consider now the behaviour of Cg{y) near the boundary y = 0. This 
is easier to analyze in the original variable r. In particular. 



Note that when C is positive, c^(r) > 1 which implies (through WKB quantization) the 
existence of the metastable states whose group velocity is greater than one [3] . Hence, for 
values of A such that C > 0, the boundary theory violates causality. That is, C should 
remain negative for the dual field theory to be consistent. The values of the Gauss-Bonnet 
parameter A for which causality is preserved are determined from the solutions of the 
inequality C < 0. To be precise. 



2 - 8A - - 16 

Note that although we have analyzed only the leading behaviour of Cg close to the boundary 
our results are exact since given by 



cl = a'{l - 4X)X{rh S7/ TeT (3-13) 



(1 - 4A) + 7A^ 
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(l-4A) + A-i (1-4A)+4A 



is a monotonically increasing function of r for all A < ^. X{r) is defined in (2.7). 

This completes the discussion of this section. The gravity analysis imposes an upper 
bound (3.12) on the Gauss-Bonnet parameter A. As we will see below, the AdS/CFT 
correspondence relates A to the coefficients of the stress-energy tensor three point function 
of the boundary CFT. Therefore, the bound on A can be translated into constraints on 
these coefficients. In the next section we will consider how similar constraints arise in field 
theory. 
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4. Energy flux one point functions and positivity of energy bounds 

Consider the integrated energy flux per unit angle measured through a very large sphere 
of radius r 

S{n) = lim r'^-^ [ dtn'T^{t,rn') (4.1) 

r— >oo J 

where n* denotes a unit vector in R^~^, the space where the field theory lives. This 
unit vector specifies the position on S^~^ where energy measurements may take place. 
Integrating over all angles yields the total energy flux at large distances. 

An interesting object to consider [1] is the energy flux one point function [30]. It is 
defined as the expectation value of the energy fliix operator (4.1) on states created by local 
operators Oq 



The simplest case to examine is when the external states Oq\0) are produced by operators 
with energy qo = q and zero momentum ,i.e., q'^ = {q, 0, 0, 0) 

O^.j^W- (4.3) 
A state with generic four momentum q^^ can be obtained by performing a simple boost on 

Oq of (4.3). 

Here we will be interested in six dimensional conformal field theories. In particular, 
we will analyze the energy fiiix one point function on states produced by the stress-energy 
tensor operator 

Oq = e,jT,j{q) (4.4) 

where is a symmetric, traceless polarization tensor with indices purely in the spatial 
directions. In complete analogy with [1] 0(5) rotational symmetry allows us to express 
(£^(n)) as 



(4.5) 

n 2 

where is the volume of the unit four-sphere Q4 = Hence, the energy fiux one point 
function is fixed by symmetry up to two coefficients, t2 and t4. 

This result is in agreement with expectations from conformal invariance. As explained 
in [31], in any d-dimensional CFT the three point function of the stress-energy tensor can 



be expressed in terms of three independent coefficients. They are denoted by a, b, c in eqs. 
(3.15) to (3.21) of [31]. On the other hand, the two point function of the stress-energy 
tensor depends on a unique parameter Ct which is related to a, b, c though 



_ 27rl {d-2){d + 3)a-2b-{d+l)c 

^ r[f] d{d + 2) ^^-^^ 

It is then convenient to "change basis" and express the three point function in terms of 
Ct and any other two hnear combinations of a, 6, c. Given that the energy flux one point 
function is actuaUy the ratio between a three- and a two- point function, it should be 
completely determined up to two independent parameters, i.e., the ratios of the two linear 
combinations of a, 6, c with Ct- 

To obtain the numbers {—\^ —^) which appear in (4.5), we require that the integral 
of the energy flux one point function over the four dimensional sphere yields the total 
energy q = Qq. To see this one can use rotational invariance to set all components of the 
polarization tensor to zero except for en = —€22- Then 

%^ = l(n? + na (4.7) 

integrated over the four sphere and divided by its volume yields |. The constant — ^ in 
the last term of (4.5) is obtained in an identical manner. 

Much like in [1], positivity of the energy fliix implies constraints on the CFT param- 
eters t2 and ^4 

l-i*.-|«4>0 

1 2 \ 4 , 

To obtain these inequalities we use rotational symmetry to set n — X5. Then e^- can be 
separated into a tensor, vector and scalar components with respect to rotations mx^ . . . x^. 
The tensor component has £5^ = e^s = €55 = and yields the flrst line in (4.8). The vector 
and scalar components give rise to the second and third line respectively. Note that each 
of the constraints (4.8) is saturated in a free fleld theory without antisymmetric tensor 
fields, fermions or scalars respectively. This is similar to the situation in four dimensions. 

It is possible to calculate the energy one point function explicitly and thus derive the 
precise relations between t4 and the coefficients which appear in the stress-energy tensor 
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three-point function a, 6, c. The computation is outhned in Appendix A. Here we present 
the results 



to = 



14 220a + 54b - 39c 



3 

t4 = -28 



36a - 26 - 7c 
65a + 246 - 14c 
36a - 26 - 7c 



t2 = -140 



2na - rif 



90na + 20?!/ + ns 

35 QUa — SUf + Us 

2 90na + 20n/ + Us 



(4.9) 



The expressions for free theories can be obtained with the help of [32]. First write a, 6, c 
of [31] in terms of A, B, C of [33] 



A 

8' 



b = 



B-2A 
8 ' 



C 

'=2 



Then use [32] to relate A,B,C with na,nf,ns ,i.e. 
two-tensors, Dirac fermions and real scalars^. 



(4.10) 

the number of free antisymmetric 



A= — 



■53' 



g3 

3 



■ 63 62 

'^S^"'* T^"'-^ 4y 6no 

■4^62 62^ ^6^ 
__n, + -8n/ + 2-6na 



(4.11) 



Note that for the supersymmetric (2, 0) multiplet with an anti-selfdual two form, 
five scalars and one Dirac fermion, ^4 vanishes identically. In general, superconformal 
Ward identities result in an additional linear constraint on the three point functions of the 
stress energy tensor which reduces the number of independent parameters to two. (See 
[34] where the explicit form of the constraint is worked out in four dimensional CFT.) 
While the precise form of the constraint is not known in six dimensions, it can be easily 
determined. Recall that supersymmetry in six dimensions implies 677-0 + — 8n/ = 
(this relation is satisfied by both scalar multiplet and (2,0) multiplet). Using (4.11) we 
can express this constraint in terms of A^ B, C as nA + 2AB — 56C = 0. This fixes the form 
of the constraint in six-dimensional superconformal theories. Note that supersymmetry in 
six dimensions implies t4 = 0, just like in four dimensions [1]. 

With ^4 = the inequalities in (4.8) reduce to 



-- < ^2 < 5 
3 - 



(4.12) 



thus constraining the domain of ^2 in any supersymmetric CFT. Even when there is 
no supersymmetry, explicit bounds on ^2,^4 can be derived from (4.8). In fact, in the 
space of t2 and ^4 the solutions of (4.8) lie within a triangle defined by the vertex points 
(-f,|),(0,f)and (7,-7). 



Note a factor of 2 missing from the last term of the last line in eq.(3.22) of [32] for d = 6. 
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5. Weyl anomaly and t2 and t4 from Gauss-Bonnet gravity 

To compare the bounds (4.8) with the causahty constraint (3.12) we need to compute t2 and 
t4. More specifically, we need to determine the coefficients a,b,c in a CFT hypothetically 
dual to GB gravity. In principle, one can compute the three-point functions of the stress- 
energy tensor directly by computing the scattering of three gravitons in the bulk of Anti 
de Sitter space. However this route is technically more challenging than the one we take 
below. Instead, we use the relation between a, 6, c and the coefficients of the Weyl anomaly. 
In a six-dimensional CFT the latter contains terms of three possible types 

3 

Aw = Ee + J2 + Vi J* (5.1) 

i 

where is the Euler density in six dimensions, = 1,...3 are three independent 
conformal invariants composed out of the Weyl tensor and its derivatives, and the last 
term is a total derivative of a covariant expression. We will use Eq and Ij in the form 
quoted in [35]. 

As explained in [35] , the second term in Aw (B-type anomaly) comes from the effective 
action which contributes to the three-point functions of the stress-energy tensor. Hence, 
there is a linear relation between the coefficients hi in (5.1) and the values of a,b,c. We 
explain how to obtain this relation below. To determine bi one needs to compute the Weyl 
anomaly in GB gravity. The procedure for computing the Weyl anomaly of a GFT in a 
holographically dual theory has been introduced in [36]. Consider a d-dimensional CFT 
formulated on a space with Euclidean metric g[^\ (We will be interested in the specific 
case d = 6.) Under a small Weyl transformation the metric changes as dg^^^ — 2dag^j\ 
The CFT action is not invariant, but rather picks up an anomalous term, 

SW[g^f] = J d'^xy/detg(^)AwSa (5.2) 

To compute the anomaly in GB gravity consider the following ansatz for the metric [36] 

where 

9io = 9^ + P9\f + + ^^P') (5-4) 

is an expansion in powers of the radial coordinate p. One can now solve the equations 
of motions of GB gravity order by order in the p expansion and determine g^j li = 1, . . . 
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in terms of g^j . One should think of gl^ as the metric which sources the stress-energy 
tensor of the boundary CFT. From the point of view of the AdS/CFT correspondence, 
specifying the source completely determines the solution (5.3). 

To compute the anomaly, one needs to substitute the resulting expansion (5.4) back 
into the lagrangian (more precisely, into y/detgC), and extract the coefficient of the 1/p 
term. This is because when integrated over p with the UV cutoff at p = e, this term gives 
rise to a loge term in the six-dimensional effective action. This term is not removed by 
local counterterms and gives rise to the anomaly (5.2) under Weyl transformation. Clearly, 
the leading term in the expansion of the lagrangian comes from y^det^ ~ l/p"^- One may 
therefore be surprised that it is suflBcient to find gij up to next-to-next-to leading order since 
an 0{p^) term in (5.4) may naively contribute to the 0{p~^) term in ^/detg C However 
one can check that this term [which is of the type g^^^^^ g^j''] contains a multiplicative factor 
which vanishes when the solution for the AdS radius (2.9) is substituted. This has been 
observed previously in the context of gravity with terms in four dimensions [37]. 

In principle, the procedure outlined above can be performed analytically to obtain 
the anomaly Aw in the form (5.1). However we opted to use the computer. Consider the 
boundary metric of the form 

6 

gijdx'dx^ = /(a;^ x^) [{dx^f + {dx^] + J^^dx'f (5.5) 

One can check that for this metric £'6 = and = 1,...3 are linearly independent 
combinations of terms which contain six derivatives of f{x^,x^) distributed in various 
ways. An example of such a term would be 

[f('^''\x^x^)]'f^'''\x^x')/[fix^x')f (5.6) 

where f^P'^^ = [d / dx'^Y{d / dx^Y f{x^ ^x'^). One can now use Mathematica to solve the 
equations of motion order by order in p. The leading non-trivial term relates the value of 
the AdS radius with the cosmological constant. The next to leading term in the equations 
of motion determines g^^\ The non- vanishing components are g'^^ — g^ , g^^ , g^^ , 
934 — 943 955 ~ 9m- '^^^ explicit expressions can be obtained by using the following 
formula 
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This is the result in Einstein-Hilbert gravity [36]. It is not modified by the inclusion of 
the finite Gauss-Bonnet term. The number of linearly independent algebraic equations at 
this order is equal to the number of the nontrivial components of glj\ At the same order, 
there are more equations which contain derivatives of g^^^ with respect to x^, x^. However 
upon substitution of (5.7) these equations are identically satisfied, which provides a good 
consistency check. This story repeats itself at the next order as well. However now A enters 

(2) 

nontrivially into the solution for g^j which is somewhat cumbersome, so we will not quote 
it here. 

The next step involves substituting (5.3) together with the solution (5.5) into the 
action (2.2) and extracting the 1/p term in the integrand. The resulting expression is too 
long to quote here, but it must be of the form J d^x\/ ^'&\.g'^^^ Aw, where Aw admits the 
representation (5.1). Both the anomaly Aw and the invariants li are long expressions 
involving terms of the type (5.6). Fortunately, the last (total derivative) term in the 
anomaly can also be represented in a convenient way. In fact, any total derivative term 
must be of the form [35] 



where Ci are certain combinations of curvature invariants which can be found in Appendix 
A of [35] , and Cj are arbitrary coefficients. Now we compute Ci for our choice of boundary 
metric (5.5) and demand that the coefficient in front of every term of the type (5.6) in 
expression 



7 




(5.8) 



3 
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(5.9) 



vanishes. This uniquely fixes hi and Ci\ the result is 




(5.10) 



70(1 + Vl -4A - 2A)(1 - 4A) 
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and 

Cl = C2 = 

C3 = 70 (3(1 + Vl - 4A) - 2A(11 + 8VI - 4A) - 20A)) 

C4 = 84 (^-3(1 + Vl - 4A) + 2A(11 + 8VI - 4A - 20A)) 

140 / / y \ (5-11) 

C5 = — (^-9(1 + Vl - 4A) + A(57 + 39Vl - 4A - 76A) j 

C6 = 14 (^-3(1 + Vl - 4A) + 2A(11 + 8VI - 4A - 20A)) 

280 / \ 
C7 = — (9{1 + Vl - 4A) - 2A(30 + 2lVl - 4A - 40A) j 

where we neglected an overall factor common to all 6j's and Cj's since only the ratios of q 
will enter the expressions for ^2 and ^4 which we are after. One can check that for A = 
the expressions (5.10) and (5.11) reduce to the numbers which appear in [35]. In this case 
the values of bi are consistent with the anomaly of a free (2,0) multiplet in six dimensions. 
The coeflBcients in front of the total derivative terms are scheme-dependent and therefore 
should not be compared. 

The coefficients bi in (5.10) are related linearly to the parameters that determine the 
two- and three-point functions of the stress-energy tensor, A,B,C. To determine this 
relation we use the free field results for the Weyl anomaly [35] : 

/28 896 8008 \ , /5 _ 2378 \, -.0^ 
Aw=- ( y"'«+~§~'^/^ — 3~""J 1 3"s-32n/ J -'2 + (2ns+40n/+180na) Is 

(5.12) 

where ns,nf,na are the same as those in eqs. (4.9) and (4.11). In (5.12) we omitted 
the A-type anomaly term since its coefficient is related to the four-point function of the 
stress-energy tensor in d — 6. We have also omitted the total derivative term in (5.12). 
Using (5.12) together with the free field results for A,B,C, (4.11), we arrive at 



A = 



864 
72 



(^3(1 + Vl - 4A) - A(25 + 19Vl -4A - 52A)) 



B = — (^181(1 + Vl - 4A) - A(1275 + 913V1 - 4A - 2204A)) (5.13) 
C = ^ (^59(1 + Vl - 4A) - A(425 + 307V1 - 4A - 756A)) 

Finally, using (4.10) and (4.9) we determine the values of a, 6, c, and ^2 and t^: 
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As explained in the previous section, vanishing is a necessary feature of any supercon- 
formal field theory. 

We are now in position to substitute the values of t2 and in (5.14) into the constraints 
(4.8) and find out what the constraints are in terms of A. The result is 

-^<A<1 (5.15) 

Note that the upper bound on A precisely coincides with the upper bound (3.12) obtained 
in Section 2 by demanding causality of the boundary theory. The lower bound in (5.15) 
presumably can be obtained from considering excitations with different polarization, just 
as it has been in the four-dimensional setup [4-5] . 



6. Discussion 

In this paper we considered Gauss-Bonnet gravity in an AdSf background. This theory has 
exact black hole solutions for values of the Gauss-Bonnet coupling A smaller than 1/4. We 
studied small fluctuations around these backgrounds and showed that causality imposes 
an upper bound (3.12) on the value of A. We expect a lower bound to follow from studying 
gravitons of different helicity. To compare the causality bound with the positivity of energy 
bounds we computed t2 and in a six-dimensional CFT in terms of the constants a, b, c 
which specify the three-point functions of the stress-energy tensor. We then computed 
the holographic Weyl anomaly, and found that t4 = 0. The resulting bounds on t2 are 
translated into bounds on A (5.15). Note that the upper bound in (5.15) is precisely the 
same as the causality bound (3.12). 

We found that the results in six dimensions are very similar to those in four dimensions. 
The fact that ^4 = is related to the truncation of the gravity action at C(i?^). It would be 
interesting to include R'^ terms and see if the nonsupersymmetric bounds can be addressed 
in this situation. A natural generalization of the Gauss-Bonnet theory to 0{R^) is third 
order Lovelock gravity. Its equations of motion again contain only terms with at most 
second derivatives acting on the metric. There are two independent coefficients which 
multiply the Gauss-Bonnet and the third order Lovelock term in the gravitational action. 
Unfortunately, at least naively, the third order Lovelock term does not contribute to the 
three-point functions of the stress-energy tensor, and hence would not affect the value of 
t4. However the situation needs to be analyzed more carefully [38]. 
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One may view our result as an additional argument in favor of the robustness of the 
correspondence between the positivity of energy and causality conditions. The understand- 
ing of this correspondence is at present somewhat incomplete although it was argued in 
[5] that for Gauss-Bonnet gravity the equation for a graviton propagating in a shock wave 
background receives contributions only from the energy flux one-point function. In fact, 
the form of eq. (3.12) suggests that the leading small temperature behaviour of the thermal 
two-point function might be determined by the three-point function at zero temperature 
alone; and it is this leading term that determines causality of the boundary theory. 

Finally, let us consider the implications of our result on the viscosity to entropy ratio. 
For a CFT hypothetically dual to d -|- 1 dimensional Gauss-Bonnet gravity, this ratio has 
been computed in [2] with the result^ 

^ = ^(i-2Aa) ^« ^ = ^(1-3A) (6.1) 



47r V d-2 J s 47r 

Combining this with the bound on A leads to 



1 7 

> — — (6.2) 



Note that the viscosity to entropy ratio is bounded from below by a number smaller than 
that in four dimensions [3]. Hence, the lower bound on viscosity in Gauss-Bonnet gravity 
depends on the dimensionality of the corresponding field theory. More precisely, the bound 
decreases as one goes from d — A to d = Q. 

This curious fact lead us to examine the restriction causality imposes on rj/s for all 
d + 1-dimensional Gauss-Bonnet theories with d < 10. The analysis is similar to the one 
done in Section 3. We find that the viscosity to entropy bound attains the smallest value 
when d = 8: 

V 
s 

It is interesting to note that the correspondence between causality and positivity of energy 
also leads to an upper bound on viscosity in the Gaus-Bonnet theories. Perhaps further 
understanding of the bounds of the type (6.3) in generalized theories of gravity may shed 
some light on the existence of a universal viscosity bound. 



1 219 , , 

> 6.3) 

d=8 ~ 47r 529 ^ ' 



Here d> 4 otherwise the Gauss-Bonnet term identically vanishes. 



15 



Acknowledgements: We thank F. Bastianelli, D. Kutasov, K. Papadodimas, L. Rastelli, 
K. Skenderis, M. Taylor, A. Tseytlin and E. Verlinde for useful discussions and correspon- 
dence. M.K. and A. P. thank the Aspen Center for Physics, the Centro de Ciencias de 
Benasque and the Simons Center for Geometry and Physics where part of this work has 
been completed. A. P. thanks the University of Amsterdam and the University of Chicago 
for hospitality. This work was partly supported by NWO Spinoza grant and the FOM 
program String Theory and Quantum Gravity. 



Appendix A. 

To determine t2 and t^ with respect to the coeflBcients which appear in the three- point 
function of the stress-energy tensor, an explicit computation of the energy fliix is necessary. 

Let us start with a careful consideration of eq. (4.5). Without loss of generality we can 
use rotational symmetry to set the detector along the X5 direction. It is then convenient to 
define new coordinates — t±X5 and express the energy fiux measured at large distances 
as ^ 

S= lim ""^ ) [ dx-T__{x+,x-) (A.l) 

a;+^oo \ 2 / J 

Since our main objective is to obtain ^2,^4 it is sufficient to extract the energy correlation 
function for two specific choices of the polarization tensor as long as they yield two inde- 
pendent linear combinations of ^2,^4- The simplest cases to consider are Sij = for all i,j 
except for 612 = €21 and Cij = for all i,j except for €15 = €51. With these two choices, 
the following linear combinations of ^2 , ^4 can be computed 

12 3 2 
1 _ _t 1 1 + _t 1 (A.2) 

5 35 ' 10 35 ^ ^ 

In what follows we will analyze the former case in detail. The latter can be treated in an 
almost identical manner. 

In practice, we need to separately consider the numerator and denominator of eq. 
(4.5). That is, we should compute the Fourier transform of the two- point function 

/2(go) = J d^xe'^^^' (Ti2(x)Ti2(0)) (A.3) 
as well as the three- point function. 

/3(go) ^ j d^xe^'^-' Juv^ (^^T^) / m2(x)T__(xi)Ti2(0)) (A.4) 
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As a warm up consider first the case of the two point function. Its form is fixed by 
conformal invariance and according to [31] can be expressed as 



(Ti2(x)Ti2(0)) = 



X 



{x^y 



(A.5) 



Here a; is a six vector parametrized as a; = {x~^ ,x^ ,1/1,112,1/3,1/4) and Ct satisfies eq. (4.6) 
with d = 6. To Fourier transform the above expression recall that the operators T12 in 
(A. 3) are ordered as written. This implies the ie prescription t ^ t — ie which in light 
cone coordinates is replaced by a;^ — > a;^ — ie. Integrating over j/i, j/2 using a spherical 
parametrization results in 

^ 2 ■ it- 28 ^^^^ ~ ~ ^""^ ^^-^^ 

where we substituted Ct in terms of the coefficients a, b, c which determine the three point 
function and denoted by I2 the integral 

We are now ready to move on to the calculation of (A.4). The starting point is once 
more the result of [31] where the form of the three- point function of the stress-energy 
tensor is determined by conformal invariance up to three numbers a, b, c. Adapting eq. 
(3.15) of [31] to the case of interest and taking the limit lim^+^^ yields 



n + 



with 

h = (a;-)' { [(48a + 86 - 7c) {y\ + y|) + c {yl + y^) + (256a + 966 - 54c) ylyl^- 

+ (48a + 86 - 6c) (y| + yf) (y^ + y|) + 2cy2y2] + ^^.9) 
+a;-a;+ [- (48a + 86 - 6c) {y\ + y^) - 2c (yf + y|)] + c (a;+)' (a;-)'} 

Note that the ie prescription here is such that the operator to the left of another acquires 
a more negative imaginary part in the time direction. When integrating over a;j~ we have 
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the option of a contour closing on either the upper or the lower Xi plane thus including 
only one of the two poles in (A. 8). This results in 

/ dx^ hm (^i^)\T^,{x,)T..{x)T,,m = ^7 ^W^^TT^ (^-10) 

J xj^oo V ^ / o {x~ — 2ie) (x^) 

Integrating now over the transverse coordinates 2/1,2/2 yields 

/3(go) = -g- (28a + 6b- 3c) h (A.ll) 

where 

Let us gather the results from eq.(A.6), (A. 7), (A.ll) and (A. 12) to form the following ratio 

^/3(go) 7 28a + 66 -3c 



Qo /2(?o) 3 36a - 26 - 7c 
As previously explained eq.(4.5) combined with (A. 13) leads to 



(A.13) 



. 1 2 7 28a + 66 -3c ^ 

^-5'- -35'^ = -3 36a -26 -7c (^'"^ 

In a similar but slightly more complicated manner it is possible to compute the other linear 
combination of t2, ^4 in (A. 2) 

3 2 ^„ 16a + 46 -3c , 

1 + —t2 14 = 28 A.15) 

10 35 36a - 26 - 7c ^ ^ 

Solving then for ^2,^4 reproduces eq. (4.9). 
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